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ABSTRACT 

Binary black-hole systems with spins aligned with the orbital angular momentum are of special 
interest as they may be the preferred end-state of the inspiral of generic supermassive binary black-hole 
systems. In view of this, we have computed the inspiral and merger of a large set of binary systems of 
equal-mass black holes with spins aligned with the orbital angular momentum but otherwise arbitrary. 
By least-square fitting the results of these simulations we have constructed two "spin diagrams" which 
provide straightforward information about the recoil velocity Vkick and the final black- hole spin a nn 
in terms of the dimensionless spins a± and ai of the two initial black holes. Overall they suggest a 
maximum recoil velocity of |i>kick| = 441.94 km/s, and minimum and maximum final spins agn = 0.3471 
and a nn = 0.9591, respectively. 

Subject headings: black hole physics - relativity - gravitational waves - stars: statistics 



1. INTRODUCTION 

A number of recent developments in numerical relativ- 
ity have allowed for stable evolution of binary black holes 
and opened the door to extended and systematic studies 
of these systems. Of particular interest to astrophysics 
are the calculations of the recoil velocity and of the spin 
of the final black hole produced by the merger. It is well 
known that a binary with unequal masses or spins will 
radiate gravitational energy asymmetrically. This results 
in an uneven flux of momentum, providing a net linear 
velocity to the final black hole. The knowledge of both 
the "kick" velocity and of the final spin could have a di- 
rect impact on studies of the evolution of supermassive 
black holes and on statistical studies on the dynamics of 
compact objects in dense stellar systems. 

Over the past year, a number of simulations have been 
carried out to determine the recoil velocities for a variety 
of binary black-hole systems. Non-spinning but unequal- 
mass binaries were the first systems to be studied and 
several works have now provided an acc urate mapping of 
the unequal-mass space of parameters ([Herrmann et aLI 
I2007at iBaker et al.ll2006bl; IGonzalez et al.ll2007b[) . More 
recently, the recoils from binaries with spinning black 
holes have also been considered by investigating equal- 
mass binaries in which the spins of the black holes 
are e ither aligned with the orbital angular mo men- 
tum ([Herrmann et alJ I2007U iKoppitz et all |2007[ ). or 
not. In the first case, a systematic investigation has 
shown that the largest recoil pos sible from such sys - 
tems is of the order of 450 km/s ([Pollnev et al.ll2007f ). 
In the second case, instead, specific configurations with 
spins orthogonal to the orbital one ha ve been shown to 
lead t o recoils as high as 250 km/s (jCampanelli et all 
I2007al IGonzalez etatl l2007aj) . suggesting a maximum 
kick o f about 4000 km/s for m aximally-spinning black 
holes (fCampa nelli et aLll2007bf ). Recoil velocities of this 
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magnitude could lead to the ejection of massive black 
holes from the hosting galaxies, with important conse- 
quences on their cosmological evolution. 

Here, we e xt end the analysis carried out 
in iPollnev et al.l (|2007f ) of binary black hole sys- 
tems with equal-mass and spins aligned with the orbital 
one. Our interest in this type of binaries stems from 
the fact that systems of this type may represent a 
preferred end-state of the binary evolution. Post- 
Newtonian studies have shown that in vacuum the 
gravitational spin-orbit coupling has a tendency to align 
the s pins when they ar e initially close to the orbital 
one ([Schnittman 2004). Furthermore, if the binary 
evolves in a disc, as expected for supermassive black 
holes, the matter can exert a tor que tending to align 
the spins (Bo gdanovic et al.ll2007l ). Finally, a recoiling 
supermassive black hole could retain the inner part of 
its accretion disc and thus the fuel for a continuing QSO 
phase lasting milli ons of year s as it moves away from the 
galactic nucleus (jLoeb Il2007|) . Yet, the analysis of QSOs 
from the Sloan Digital Sky Survey shows no evidence 
for black holes carrying an ac cretion disc and hen ce for 
very large recoiling velocities ([Bonning et al.ll2007f ). 

2. NUMERICAL SETUP AND INITIAL DATA 

The numerical simulations have been carried out 
using the CCATIE code, a three-dimensional finitc- 
differenc es code usin g the Cactus Computational 
Toolkit (ICactus II2007D and Carpet mesh refinement in- 
frastructure (|Schnetter et ai1l2007| ). The main features 
of the code have been recently reviewed in IPollnev et al.l 
(2007), where the code has been employe d using the 
so-called "moving-punctures" technique (jBaker et al.l 
l2006at ICampanelli et all 12006a) . The initial data con- 
sists of five sequences with constant orbital angular mo- 
mentum, which is however different from sequence to se- 
quence. In the r and ra-sequences, the initial spin of 
one of the black holes S2 is held fixed along the z-axis 
and the spin of the other black hole is varied so that 
the spin ratio a 1/02 takes the values between —1 and 
+ 1, with a, = Si/M/. In the i-sequence, instead, the 
spin with a negative z-component is held fixed, while 
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TABLE 1 

Binary sequences for which numerical simulations have been carried out, with different columns referring to the 
puncture initial location ±x/m, the linear momenta ±p/m, the mass parameters itli/m, the dimensionless spins oj, the 

normalized adm mass m adm = m adm /m measured at infinity, and the normalized adm angular momentum 

J ADM = J ADM /M 2 . Finally, the last six columns contain the numerical and fitted values for |«kickl ( IN km/s), aa n and the 

CORRESPONDING ERRORS. 



±x/M ±p/M m 1 /M m 2 /M 



o 2 



M 5 



J. 



b fit I 

1 kigk 1 



err. (%) a_ 



tin 



fit 

iiL. 



a 



err. (%) 



rO 
rl 
r2 
r3 
rl 
r5 
;-6 
rl 
r8 



3.0205 
3.1264 
3.2198 
3.3190 
3.4100 
3.5063 
3.5988 
3.6841 
3.7705 



0.1366 
0.1319 
0.1281 
0.1243 
0.1210 
0.1176 
0.1146 
0.1120 
0.1094 



0.4011 
0.4380 
0.4615 
0.4749 
0.4796 
0.4761 
0.4638 
0.4412 
0.4052 



0.4009 
0.4016 
0.4022 
0.4028 
0.4034 
0.4040 
0.4044 
0.4048 
0.4052 



-0.584 
-0.438 
-0.292 
-0.146 
0.000 
0.146 
0.292 
0.438 
0.584 



0.584 
0.584 
0.584 
0.584 
0.584 
0.584 
0.584 
0.584 
0.584 



0.9856 
0.9855 
0.9856 
0.9857 
0.9859 
0.9862 
0.9864 
0.9867 
0.9872 



0.825 
0.861 
0.898 
0.935 
0.971 
1.007 
1.044 
1.081 
1.117 



261.75 
221.38 
186.18 
144.02 
106.11 
81.42 
45.90 
20.59 
0.00 



258.09 
219.04 
181.93 
146.75 
113.52 
82.23 
52.88 
25.47 
0.00 



1.40 
1.06 
2.28 
1.90 
6.98 
1.00 
15.21 
23.70 
0.00 



0.6891 
0.7109 
0.7314 
0.7516 
0.7740 
0.7948 
0.8150 
0.8364 
0.8550 



0.6883 
0.7105 
0.7322 
0.7536 
0.7747 
0.7953 
0.8156 
0.8355 
0.855 



0.12 
0.06 
0.11 
0.27 
0.08 
0.06 
0.07 
0.11 
0.00 



raO 

ral 

ra2 

ra3 

ra4 

ra6 

ra8 

ralO 

ral2 



2.9654 
3.0046 
3.0438 
3.0816 
3.1215 
3.1988 
3.2705 
3.3434 
3.4120 



0.1391 
0.1373 
0.1355 
0.1339 
0.1321 
0.1290 
0.1261 
0.1234 
0.1209 



0.4585 
0.4645 
0.4692 
0.4730 
0.4757 
0.4782 
0.4768 
0.4714 
0.4617 



0.4584 
0.4587 
0.4591 
0.4594 
0.4597 
0.4602 
0.4608 
0.4612 
0.4617 



-0.300 
-0.250 
-0.200 
-0.150 
-0.100 
0.000 
0.100 
0.200 
0.300 



0.300 
0.300 
0.300 
0.300 
0.300 
0.300 
0.300 
0.300 
0.300 



0.9845 
0.9846 
0.9847 
0.9848 
0.9849 
0.9850 
0.9852 
0.9853 
0.9855 



0.8250 
0.8376 
0.8499 
0.8628 
0.8747 
0.9003 
0.9248 
0.9502 
0.9750 



131.34 
118.10 
106.33 
94.98 
84.74 
63.43 
41.29 
19.11 
0.00 



132.58 
120.28 
108.21 
96.36 
84.75 
62.19 
40.55 
19.82 
0.00 



0.95 
1.85 
1.77 
1.46 
0.01 
1.95 
1.79 
3.72 
0.00 



0.6894 
0.6971 
0.7047 
0.7120 
0.7192 
0.7331 
0.7471 
0.7618 
0.7772 



0.6883 
0.6959 
0.7035 
0.7111 
0.7185 
0.7334 
0.7481 
0.7626 
0.7769 



0.16 
0.17 
0.17 
0.13 
0.09 
0.04 
0.13 
0.11 
0.03 



sO 
si 
s2 
s3 
s4 
s5 



2.9447 
3.1106 
3.2718 
3.4098 
3.5521 
3.6721 
3.7896 



0.1401 
0.1326 
0.1261 
0.1210 
0.1161 
0.1123 
0.1088 



0.4761 
0.4756 
0.4709 
0.4617 
0.4476 
0.4276 
0.4002 



0.4761 
0.4756 
0.4709 
0.4617 
0.4476 
0.4276 
0.4002 



0.000 
0.100 
0.200 
0.300 
0.400 
0.500 
0.600 



0.000 
0.100 
0.200 
0.300 
0.400 
0.500 
0.600 



0.9844 
0.9848 
0.9851 
0.9855 
0.9859 
0.9865 
0.9874 



0.8251 
0.8749 
0.9251 
0.9751 
1.0250 
1.0748 
1.1246 



0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



0.6892 
0.7192 
0.7471 
0.7772 
0.8077 
0.8340 
0.8583 



0.6883 
0.7185 
0.7481 
0.7769 
0.8051 
0.8325 
0.8592 



0.13 
0.09 
0.13 
0.03 
0.33 
0.18 
0.11 



to 
fl 

£2 
f3 
U 
t5 
tt) 

a 
ts 



4.1910 
4.0812 
3.9767 
3.8632 
3.7387 
3.6102 
3.4765 
3.3391 
3.1712 



0.1074 
0.1103 
0.1131 
0.1165 
0.1204 
0.1246 
0.1294 
0.1348 
0.1419 



0.4066 
0.4062 
0.4057 
0.4053 
0.4047 
0.4041 
0.4033 
0.4025 
0.4015 



0.4064 
0.4426 
0.4652 
0.4775 
0.4810 
0.4761 
0.4625 
0.4387 
0.4015 



-0.584 
-0.584 
-0.584 
-0.584 
-0.584 
-0.584 
-0.584 
-0.584 
-0.584 



0.584 
0.438 
0.292 
0.146 
0.000 
-0.146 
-0.292 
-0.438 
-0.584 



0.9889 
0.9884 
0.9881 
0.9879 
0.9878 
0.9876 
0.9874 
0.9873 
0.9875 



0.9002 
0.8638 
0.8265 
0.7906 
0.7543 
0.7172 
0.6807 
0.6447 
0.6080 



259.49 
238.37 
200.25 
174.58 
142.62 
106.36 
71.35 
35.36 
0.00 



258.09 
232.62 
205.21 
175.86 
144.57 
111.34 
76.17 
39.05 
0.00 



0.54 
2.41 
2.48 
0.73 
1.37 
4.68 
6.75 
10.45 
0.00 



0.6868 
0.6640 
0.6400 
0.6180 
0.5965 
0.5738 
0.5493 
0.5233 
0.4955 



0.6883 
0.6658 
0.6429 
0.6196 
0.5959 
0.5719 
0.5475 
0.5227 
0.4976 



0.22 
0.27 
0.45 
0.26 
0.09 
0.33 
0.32 
0.11 
0.42 



til 
tt2 
u3 
tt 1 



2.9500 
2.9800 
3.0500 
3.1500 



0.1398 
0.1384 
0.1355 
0.1310 



0.4683 
0.4436 
0.3951 
0.2968 



0.4685 
0.4438 
0.3953 
0.2970 



-0.200 
-0.400 
-0.600 
-0.800 



0.200 
0.400 
0.600 
0.800 



0.9845 
0.9846 
0.9847 
0.9850 



0.8248 
0.8249 
0.8266 
0.8253 



87.34 
175.39 
266.39 
356.87 



88.39 
176.78 
265.16 
353.55 



1.20 
0.79 
0.46 
0.93 



0.6893 
0.6895 
0.6884 
0.6884 



0.6883 
0.6883 
0.6883 
0.6883 



0.15 
0.17 
0.01 
0.01 



in the s and u-sequences a\jai = 1 and —1, respec- 
tively. In all cases, the masses are Mi = M/2 = 1/2. For 
the orbital initial data parameters we use the effective- 
potential method, which allows one to choose the initial 
data parameters such that the resulting physical parame- 
ters {e.g., masses and spins) describe a binary black- hole 
system on a quasi-circular orbit. The free parameters 
are: the coordinate locations the mass parameters 
rrii, the linear momenta Pi, and the spins S^. Quasi- 
circular orbits are then selected by setting pi = — P2 to 
be orthogonal to C2 — Ci, so that L = Ci x pi + C2 x P2 
is the orbital angular momentum. The initial parameters 
are collected in the left part of Table HJ while the right 
part reports the results of simulations. For all of them we 
have employed 8 levels of refinement and a minimum res- 
olution 0.024 M, which has been reduced to 0.018 M for 
binaries r5, r6. Note that our results fo r the u-sequence 
differ s lightly from those reported by iHerrmann et al.l 
(|2007bf ). probably because of our accounting of the in- 



tegration constant in |wki c k| ()Pollnev et alj |2007). 

3. SPIN DIAGRAMS AND FITS 

Clearly, the recoil velocity and the spin of the final 
black hole are among the most important pieces of infor- 
mation to be extracted from the inspiral and coalescence 
of binary black holes. For binaries with equal masses 
and aligned but otherwise arbitrary spins, this informa- 
tion depends uniquely on the dimensionless spins of the 
two black holes 01, 02 and can therefore be summarized 
in the portion of the (ai, 02) plane in which the two spins 
vary. It is therefore convenient to think in terms of "spin 
diagrams" , which summarize in a simple way all of the 
relevant information. In addition, since the labelling "1" 
and "2" is arbitrary, the line a\ — 0,2 in the spin diagram 
has important symmetries: the recoil velocity vector un- 
dergoes a 7r-rotation, i.e., «kick(ai, 0-2) = — Wkick(a2, eti) 
but |ukick(ai, 02)! = |wkick(a2, d\ ) \ , while no change is ex- 
pected for the final spin, i.e., afi n (ai, CL2) = »fm(a2, a±). 
These symmetries not only allow us to consider only one 



Spin Diagrams for bh binaries with aligned spins 
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Fig. 1. — Position in the (oi, 122) space of the five sequences 
r, ra, s, t, and u for which the inspiral and merger has been com- 
puted. 

portion of the (01, 02) space (c/. Fig. [I}, thus halving 
the computational costs (or doubling the statistical sam- 
ple), but they will also be exploited later on to improve 
our fits. The position of the five sequences within the 
(oi) 02) space is shown in Fig. Q] 

Overall, the data sample computed numerically con- 
sists of 38 values for |t>kick| and for aa n which, for sim- 
plicity, we have considered to have constant error-bars 
of 8 km/s and 0.01, whic h represent, r espec tively, the 
largest errors reported in iPollnev et al.l (|2007f ). In both 
cases we have modelled the data with generic quadratic 
functions in a\ and 02 so that, in the case of the recoil 
velocity, the fitting function is 

hdckl = c + ciai + c 2 a\ + d aia 2 + d 1 a 2 + d 2 a\\ . (1) 

Note that the fitting function on the right-hand-side of 
([1} is smooth everywhere but that its absolute value is 
not smooth along the diagonal a\ — a 2 . Using |T]) and a 
blind least-square fit of the data, we obtained the coeffi- 
cients (in km/s) 

c = 0.67 ± 1.12, d = -18.56 ±5.34, 
ci = -212.85 ±2.96, d l = 213.69 ± 3.57 , 

c 2 = 50.85 ±3.48, d 2 = -40.99 ± 4.25 , (2) 

with a reduced-^ 2 = 0.09. Clearly, the errors in the co- 
efficients can be extremely large and this is simply the 
result of small-number statistics. However, the fit can 
be improved by exploiting some knowledge about the 
physics of the process to simplify the fitting expressions. 
In particular, we can use the constraint that no recoil 
velocity should be produced for binaries having the same 
spin, i.e., that |t>kick| = for a\ = a 2l or the symme- 
try condition across the line a\ — a 2 . Enforcing both 
constraints yields 

c = , ci = -di , c 2 = -d 2 , d Q = , (3) 
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Fig. 2. — Contour plots of l^kickl as a function of the spin pa- 
rameters ai and 02. The diagram has been computed using ex- 
pressions l|4} and J5j. 

thus reducing the fitting function (JTJ) to the simpler ex- 
pression 

Kickl = |ci(ai - a 2 ) + c 2 (al - a|)| . (4) 

Performing a least-square fit using Q we then obtain 

ci = -220.97 ±0.78, c 2 = 45.52 ± 2.99 , (5) 

with a comparable reduced-x 2 = 0.14, but with error- 
bars that are much smaller on average. Because of this, 
we consider expression Q as the best description of the 
data at second-order in the spin parameters. Using (01 
and ([5]), we have built the contour plots shown in Fig.[2j 
A few remarks are worth making. Firstly, we recall 
that post-Newtonian calculations have so far derived only 
the linear contribution i n the spin to the recoil velocity 
(see iFavata et al.l (|2004[) and references therein) . How- 
ever, the size of the quadratic coefficient ([5]) is not small 
when compared to the linear one and it can lead to 
rather sizeable corrections. These are maximized when 
ai = and a 2 = ±1, or when a\ = ±1 and a 2 = 0, 
and can be as large as ~ 20%; while these corrections 
are smaller than those induced by asymmetries in the 
mass, they are instructive in pointing out the relative 
importance of spin-spin and spin-orbit effects during the 
merger and can be used as a guide in further refine- 
ments of the post-Newtonian treatments. Secondly, ex- 
pression (j4|) clearly suggests that the maximum recoil ve- 
locity should be found when the asymmetry is the largest 
and the spins are antiparallel, i.e., a\ = —a 2 . Thirdly, 
when a 2 = const., e xpression (|4l) confirm s the quadratic 
scaling proposed in IPollnev et al.l (|2007h with a smaller 
data set [cf., eq. (42) there]. Fourthly, for ai = — a 2l 
expression (|4]) is only linear and reprod uces the scal- 
ing suggested by iHerrmann et alJ (|2007bl ) . Finally, us- 
ing ([J) the maximum recoil velocity is found to be 
l^kickl = 441. 94± 1.56 km/s, in very good agreement with 
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the re sults of lHerrmann et all (|2007bT) and lPollnev et al.1 

(Hqq3). 

In the same way we have first fitted the data for dfi n , 
with a function 

dfin = Pq + Pxax + p 2 a\ + qoaia 2 + q\a 2 + q 2 a\ , (6) 

and found coefficients with very large error-bars. As a 
result, also for aa n we resort to physical considerations 
to constrain the coefficients po . . . q 2 . More specifically, 
we expect that, at least at lowest order, binaries with 
equal and opposite spins will not contribute to the final 
spin and thus behave essentially as nonspinning binaries. 
Stated differently, we assume that ag n = po for binaries 
with a\ = —a 2 . In addition, enforcing the symmetry 
condition across the line a\ — a 2 we obtain 

pi =q%, P2=q2 = qo/2, (7) 

so that the fitting function ([6]) effectively reduces to 

flfin = Po +P\{a\ + a 2 ) +p 2 (ai + a 2 ) 2 . (8) 

Performing a least-square fit using {5| we then obtain 

pa = 0.6883 ± 0.0003 , p x = 0.1530 ± 0.0004 , 

p 2 = -0.0088 ±0.0005, (9) 

with a reduced-x 2 = 0.02. 

It should be noted that the coefficient of the quadratic 
term in @ is much smaller then the linear one and with 
much larger error-bars. Given the small statistics it is 
hard to assess whether a quadratic dependence is nec- 
essary or if a linear one is the correct one (however, 
see also the comment below on a possible interpreta- 
tion of expression ([8])). In view of this, we have re- 
peated the least-square fit of the data enforcing the con- 
ditions ([7]) together with p 2 = {i.e., adopting a linear 
fitting function) and obtained po = 0.6855 ± 0.0007 and 
px = 0.1518 ± 0.0012, with a worse reduced-x 2 = 0.16. 
Because the coefficients of the lowest-order terms are so 
similar, both the linear and the quadratic fits are well 
within the error-bars of the numerical simulations. Nev- 
ertheless, since a quadratic scaling yields smaller resid- 
uals, we consider it to be the best representation of the 
data and have therefore computed the contour plots in 
Fig. [3] using © and ©. 

Here too, a few remarks are worth making: Firstly, the 
fitted value for the coefficient po a grees very well with the 
values reported by several groups ((Gonzalez et al.ll2007bl ; 
iBerti et al. Il2007n when studying the inspiral of unequal- 
mass nonspinning binaries. Secondly, expression ([5]) has 
maximum values for a\ — a 2 , suggesting that the max- 
imum and minimum spins are Ofi n = 0.9591 ± 0.0022 
and a fin = 0.3471 ± 0.0224, respectively. Thirdly, the 
quadratic scaling for qfi n sub s tantial ly confirms the sug- 
gestions of ICampanelli et all (|2006bl ) but provides more 
accurate coefficients. Finally, although very simple, ex- 
pression ^ lends itself to an interesting interpretation. 
Being effectively a power series in terms of the initial 
spins of the two black holes, its zeroth-order term can 
be seen as the orbital angular momentum not radiated 
in gravitational waves and which amounts, at most, to 
~ 70% of the final spin. The first-order term, on the 
other hand, can be seen as the contribution to the final 
spin coming from the initial spins of the two black holes 
and this contribution, together with the one coming from 
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Fig. 3. — Contour plots of aa n as a function of the spin param- 
eters ai and a,2- The diagram has been computed using expres- 
sions JHJ and J5}. 

the spin-orbit coupling, amounts at most to ~ 30% of the 
final spin. Finally, the second-order term, which is natu- 
ral to expect as nonzero in this view, can then be related 
to the spin-spin coupling, with a contribution to the final 
spin which is of ~ 4% at most. 

As a side remark we also note that the monotonic 
behaviour expressed by ([9]) does not show the pres- 
ence of a local maximum of aa n c ± 0.87 for a% — 
a 2 ~ 0.34 as suggested by iDamour 1 (|2001[ ) in the ef- 
fective one-body (EOB) approximation. Because the 
latter has been shown to be in good agreement with 
numer ical-relativity simulations of nonspinning blac k 
holes (|Damour and Nagar I [20071 IDamour et al. Il2007h . 
additional simulations will be necessary to refute these 
results or to improve the EOB approximation for spin- 
ning black holes. 

Reported in the right part of Table Q] are also the fit- 
ted values for a nn and |i>kick| obtained through the fitting 
functions ([4]) and §8\), and the corresponding errors. The 
latter are of few percent for most of the cases and increase 
up to ~ 20% only for those binaries with very small kicks 
and which are intrinsically more difficult to calculate. As 
a concluding remark we note that the fitting coefficients 
computed here have been constructed using overall mod- 
erate values of the initial spin; the only exception is the 
binary uA which has the largest spin and which is never- 
theless fitted with very small errors (cf. Table [l}. In ad- 
dition, since the submission of this work, another group 
has reported results from e qual-mass binaries with spins 
as high as a x = a 2 = ±0.9 (jMarronetti et all [20071 ). Al- 
though also for these very high-spin binaries the error in 
the predicted values is of 1% at most, a larger sample 
of high-spin binaries is necessary to validate that the fit- 
ting expressions (TJ]) and © are robust also at very large 
spins. 
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4. CONCLUSIONS 

We have performed least-square fits to a large set of 
numerical-relativity data. These fits, combined with 
symmetry arguments, yield analytic expressions for the 
recoil velocity and final black hole spin resulting from the 
inspiral and merger of equal-mass black holes whose spins 
are parallel or antiparallel to the orbital angular momen- 
tum. Such configurations represent a small portion of the 
space of parameters, but may be the preferred ones if 
torques are present during the evolution. Using the ana- 
lytic expressions we have constructed two spin diagrams 
that summarize simply this information and predict a 
maximum recoil velocity of |i>kick| = 441.94 ± 1.56 km/s 
for systems with ai = —a 2 = 1 and maximum (mini- 
mum) final spin a fi „ = 0.9591 ± 0.0022 (0.3471 ± 0.0224) 
for systems with a\ = <Z2 = 1 (— 1). 
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NOTE ADDED IN PROOF. Since the publication on 
the preprint archive of this analysis, our work on the 
modelling of the final spin has progressed rapidly, yield- 
ing new results that complement and complete the 
o nes presented h e re. In particular, the work published 
in lRezzolla et al.l (|2008af) complements the analysis car- 
ried here to unequal-mass, equal-spin aligned bina ries, 
while the work reported in iRezzolla et alJ (|2008bl ) ex- 
tends it to generic binaries. 
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